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Abstract. We present a phenomenological study of charm quark photoproduction in longitudinally po- 
larized lepton-hadron collisions at next-to- leading order accuracy of QCD. Our results are based on a 
recently developed, flexible parton-level Monte Carlo program for spin-dependent heavy flavor hadropro- 
duction, which we extend to deal also with both direct and resolved photon contributions. The subsequent 
hadronization into charmed mesons is modeled in our calculations, which allows us to compare with data 
on double-spin asymmetries for D° meson production taken by the COMPASS collaboration. In general, 
next-to-leading order QCD corrections are found to be very significant and do not cancel in spin asym- 
metries. We elucidate the role of the individual hard scattering subprocesses and determine the range of 
parton momentum fractions predominantly probed for charm production at COMPASS. Theoretical un- 
certainties are estimated by varying renormalization and factorization scales and parameters controlling 
the hadronization of the charm quarks. 
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1 Motivation and Introduction 

The quest to understand the partonic structure of the nu- 
cleon spin remains to be one of the key research areas in 
Hadronic Physics even after more than 25 years of strenu- 
ous experimental and theoretical efforts. In particular, the 
polarized gluon density, defined as 



Ag(x,fi) = g+(x,fj.) - g-(x,(i) 



(1) 



is still one of the most elusive quantities associated with 
the non-perturbative partonic structure of hadrons. Here, 
9+ (.9-) denotes the probability of finding a gluon at a 
scale /i with light-cone momentum fraction x and helicity 
+ (— ) in a proton with helicity +. The total, x integrated 
gluon polarization, 



Ag(v) = / Ag(x,jj) dx 



(2) 



and a similar contribution from the sum of all quarks and 
antiquarks, enters the helicity sum rule of the nucleon 
along with the orbital angular momenta of quarks and 
gluons pp. The challenge is to precisely map Ag(x,fi) in 
a wide range of x in order to minimize extrapolation un- 
certainties in the first moment @. 

Currently, the best constraints on Ag{x, /i) are de- 
rived from global QCD analyses [HE] which treat all avail- 
able experimental probes simultaneously and consistently 
at a given order in the strong coupling a s in perturba- 
tive QCD. The availability of next-to-leading order (NLO) 



QCD corrections is essential for any meaningful, quanti- 
tative analysis of parton densities. In particular, single- 
inclusive pion [4 an d jet production [5] , measured in spin- 
dependent proton-proton collisions at BNL-RHIC [23 , have 
started to put significant limits on the amount of gluon 
polarization in the nucleon [2]. New, preliminary single 
and di-jet data from the STAR collaboration '!7\ show for 
the first time tantalizing hints for a non-zero Ag(x, //) 
8,9 . Due to the given kinematics, the current probes 
mainly constrain Ag in the medium-to-large x region, 
0.05 < x < 0.2, which is not sufficient to reliably deter- 
mine its integral © . A very significant contribution of up 
to one unit of h, i.e., twice the proton spin, can still come 
from the unexplored small x region Narrowing down 
the uncertainties on Ag(x,/i) and, at the same time, ex- 
tending the range in x continues to be the main objective 
of experimental efforts in the years to come. Eventually, 
only a future high-energy polarized electron-proton col- 
lider, such as the EIC project [TU], will finally be able to 
quantitatively address all the remaining open questions 
related to the helicity structure in the small x region [llj . 

World-data on polarized inclusive and semi-inclusive 
deep-inelastic scattering are pivotal in constraining the 
helicity quark and antiquark densities [2HB]) but due to 
lack of sufficient kinematic coverage of present fixed-target 
experiments, information on the gluon density is noto- 
riously difficult to obtain from QCD scaling violations. 
Viable probes of Ag(x,fi) at fixed-target energies com- 
prise of one- and two-hadron and open charm production 
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Fig. 1. Approximate LO extractions of Ag/g from polarized 
lepton-nucleon scattering experiments [12J131 1 14115] compared 
to results of the DSSV global analysis of helicity PDFs [2] based 
on RHIC pp data for two different scales Q. 



and have been exploited by several experimental collabo- 
rations [miT5irHlll5) . The proper theoretical description 
of these processes depends on the virtuality of the probing 
photon and is in general more involved than correspond- 
ing calculations for hadron-hadron scattering. In case of 
photoproduction, where a quasi-real photon is exchanged, 
one has to include also "resolved" contributions, where the 
photon fluctuates into a vector meson of the same quan- 
tum numbers before the hard scattering with partons in 
the proton takes place. If the virtuality Q of the photon is 
of 0(1 GeV) or higher, resolved processes are sufficiently 
suppressed but the additional momentum scale Q greatly 
complicates the calculations of phase-space and loop inte- 
grals. 

Not surprisingly, only very few calculations at NLO of 
QCD are available in case of polarized beams and targets 
[T5UT71IT8] . none of which for large virtualities Q. As a 
consequence, the available data sets on hadron and open 
charm production [T!^[T5llT4"Ill"5] have not been included 
in global QCD analyses of helicity parton densities so far. 
Experiments have analyzed their data only in terms of the 
gluon polarization, Ag(x, fi)/g{x, fi), under certain simpli- 
fying assumptions and based on leading order (LO) ma- 
trix elements. Nonetheless, the results of these exercises, 
illustrated in Fig. [TJ are in fairly good agreement with a 
NLO extraction of Ag(x, /it) including the RHIC pp data 
[2 r While determinations of Ag at hxed-target energies 
are more involved and provide less of a constraint than 
collider data, they are crucial for further testing and es- 
tablishing the assumed universality of helicity-dependent 
parton densities and hence for our understanding of the 
spin structure of the nucleon and QCD in general. On the 
one hand, lepton-nucleon scattering experiments are sen- 
sitive to different partonic hard scattering processes than 



jet or hadron production at RHIC and, on the other hand, 
the relevant momentum fractions x fall within the range 
already probed by pp data. 

In this paper, we present a comprehensive phenomcno- 
logical analysis of longitudinally polarized photoproduc- 
tion of heavy flavors at NLO of QCD. To this end, we 
extend our recently developed flexible parton-level Monte 
Carlo program for spin-dependent heavy flavor hadropro- 
duction [T2] by including all relevant subprocesses for di- 
rect [TfjlfTT] and resolved ,19,20] photon contributions. To 
facilitate the comparison to data, the hadronization of the 
produced charm quarks into charmed mesons and, as an 
additional option, their semi-leptonic decays can be mod- 
eled in our calculations based on the phenomenological 
functions used in Ref. [21] . The presented results will al- 
low one to consistently include available and upcoming 
data on spin-dependent open charm photoproduction into 
future global QCD analyses of helicity parton densities at 
NLO accuracy. 

As a first application, we examine in detail all aspects 
of open charm production in the kinematic regime accessi- 
ble with the COMPASS experiment [TS]. We will demon- 
strate that NLO corrections are very significant and do 
not cancel in double-spin asymmetries used to extract in- 
formation on the polarized gluon density Ag(x y /i). We 
elucidate the relevance of the individual direct and re- 
solved partonic subprocesses and estimate theoretical un- 
certainties by varying renormalization and factorization 
scales and parameters controlling the hadronization of the 
produced charm quarks into experimentally observed D 
mesons. Finally, we shall illustrate the sensitivity of the ex- 
isting data on the double-spin asymmetry for charm quark 
photoproduction 15,22,23,24 to Ag(x,fi) and compare 
our estimates of the relevant range of momentum fractions 
x with those obtained by the COMPASS collaboration. 

We note that our flexible Monte Carlo code is capable 
of computing any infrared safe heavy flavor photoproduc- 
tion cross section at 0(aa 2 s ) in longitudinally polarized 
lepton-nucleon collisions, including correlations of the pro- 
duced heavy quark pair. Our results complement and sig- 
nificantly extend previously existing spin-dependent NLO 
calculations of single-inclusive heavy quark yields based 
on largely analytical methods [TBlll7j . where the resolved 
contribution was neglected, any information on the par- 
tonic recoil system was lost, and most experimental cuts 
could not be implemented. The results presented here will 
be also useful for future studies of spin-dependent charm 
photoproduction at an EIC [10] . 

The outline of the paper is as follows: in Sec. II we 
briefly review some of the technical aspects of setting up a 
parton-level Monte Carlo program for heavy flavor photo- 
production in polarized lepton-nucleon collisions at NLO 
accuracy. In Sec. Ill we present a detailed phenomenolog- 
ical study of open charm quark production at COMPASS, 
including the relevance of NLO corrections, the role of 
the different hard scattering subprocesses, an estimate of 
the relevant momentum fractions x, and an assessment 
of theoretical uncertainties. We summarize our results in 
Sec. IV. 
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2 Technical Framework 

We consider heavy quark photoproduction in longitudi- 
nally polarized lepton-nucleon collisions at NLO accuracy 
of QCD by consistently including for the first time both 
direct and resolved photon contributions. In order to com- 
pute arbitrary infrared-safe observables within flexible ex- 
perimental acceptance cuts and to account for the hadroni- 
zation of the produced heavy quark pair, all phase-space 
integrations are performed numerically with appropriate 
Monte Carlo techniques. To this end, we follow closely 
the subtraction method devised and used in Refs. [251126] . 
In the following, we will only briefly outline the technical 
details relevant for the case of polarized photoproduction. 

Assuming, as usual, factorization, the inclusive cross 
section for producing a heavy quark Q in spin-dependent 
lepton-proton collisions at a center-of-mass system (c.m.s.) 
energy y/S can be written as a convolution, 

dAa Q ee i [d<rj + - d<7®_ 

= J dx a dx b Af l a (x a ,n f )AfP(x b ,n f )S 

a, b 

x dAa ab (xa,Xb,S,mQ,ki,k 2 ,(if,Hr), (3) 

where the subscripts ± in © label the helicity states of the 
lepton and proton. In analogy to Eq. ([T]), the Af a '^ denote 
the spin-dependent parton distribution functions (PDFs) 
of flavor a, b in a lepton I or proton p. The sum in ([3]) is 
over all contributing partonic processes ab — > QQc, includ- 
ing those with a direct photon, i.e., = 7, with dAa ab the 
associated polarized hard scattering cross sections. The re- 
quired spin-dependent matrix elements squared Z\|M a (,| 2 
at NLO accuracy in d — 4 — 2s dimensional regularization 
for the direct and resolved photon processes can be taken 
from [16] and [20], respectively. Corresponding unpolar- 
ized results can be found in Refs. [23112^1127] . At 0(aa 2 s ), 
parton c can be either a gluon or a light (anti-)quark pro- 
ducing an additional jet or hadron which is usually not 
observed in experiment. ki t 2 denote the four-momenta of 
the heavy quark Q and antiquark Q with mass tuq, i.e., 
fc[ 2 = m Q- Mr and fif are the renormalization and factor- 
ization scale, respectively, which are typically chosen as a 
combination of the hard scales characterizing the process. 

The "measurement function" S in Eq. ^ defines the 
observable, for instance, through a set of step functions 
implementing the experimental cuts imposed on the final- 
state particles and selecting a certain bin in a kinematical 
distribution of interest. In case of the COMPASS exper- 
iment, single-inclusive spectra of D mesons, differential 
in transverse momentum p®, are measured |15il221l2"3"ll24j . 
Thus, the cross section ^ at the heavy quark- level is not 
yet sufficient for comparing theory with experimental re- 
sults and needs to be convoluted with an additional phe- 
nomenological function D^^ Hq modeling the hadroniza- 
tion of the charm quark: 

dAo H » = dAo Q ® D Q ^ Bc > . (4) 



We will specify our choice for D®^ Hq , along with all other 
non-perturbative inputs, in Sec. [3] Even though not rele- 
vant for the phenomcnological applications considered in 
this paper, our flexible Monte Carlo program can also ac- 
count for the semi-leptonic decays of D mesons if needed. 

For the integration of the fully exclusive partonic cross 
sections dAa ab in Eq. ([3]), we generalize the framework 
used for unpolarized heavy quark production in 25,25] to 
deal with singular regions of phase-space. The gist of the 
method is to add and subtract appropriate spin-dependent 
counter terms to dAa ab which can be integrated analyti- 
cally with respect to momenta of unresolved partons. For 
all resolved photon processes, o / 7, we can adopt the 
expressions obtained in our recent calculation of polarized 
heavy quark hadroproduction at O{o 3 s ) [TJJ], since these 
processes share the same hard scattering matrix elements. 
Hence, we only need to consider processes with a direct 
photon, dA<7y b , in the following. 

In general, for a numerically efficient implementation 
of the subtraction method it is convenient 25,26,19 to 
express the three-body phase-space in terms of variables 
where soft and collinear singularities are particularly trans- 
parent. This is achieved by choosing x = (fci + fc2) 2 /s, the 
invariant mass of the QQ pair scaled by the available par- 
tonic c.m.s. energy squared, i.e., 4toq/s < x < 1. In addi- 
tion, one uses — 1 < y < 1, the cosine of the angle between 
the z-direction, aligned with the spatial direction of par- 
ton a, and fc3, the momentum of parton c, in the c.m.s. of 
the incoming partons. Soft and collinear regions of phase- 
space are then simply associated with x = 1 and y = ±1, 
respectively, and the hard scattering matrix element for 
the process 76 — > QQc can be written as 

^iM 7b |^ ^ ( / m ^;: y ' 01 2 f 2) , (5) 

s z {l — x)- ! {l — y z ) 

where Af lb is regular for x = 1 and y = ±1. The angles 
01,2 are used to parametrize the spatial orientation of fci,2 
with respect to the plane span by the other three momenta 
in the c.m.s. of the QQ pair; for further details, see [55,26, 
We note that the genuine NLO subprocess 79(g) — > 
QQq(q) can have only collinear singularities at 0(aa 2 s ). 

The partonic subprocesses contributing to the direct 
photon cross section at 0(aa 2 ), 75 — > QQg and jq (q) — > 
QQq (<?), can be decomposed as [26] 

dAo lb = dAa® + dAo^ b +) + dAa^ 

+ dAa\ s) b +dAo$ +dAa i i { ) . (6) 

Here, dAa^ and dAaty denote the 0(aa s ) Born contri- 
bution and the 0(aa 2 s ) one-loop corrections to the 75 — > 
QQ process, respectively. Analytic expressions for the vir- 
tual contributions in d dimensions, with ultraviolet diver- 
gences being subtracted at the renormalization scale /i r , 
can be found in Ref. [16]. 

In Eq. ((6]), dAa^l denotes the soft gluon emission part 
of the 7<jf scattering cross section, which can be obtained 
from the full <i-dimensional matrix elements squared in 
the limit x 1 where the phase-space integrations can 
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be performed analytically [25] . We obtain the same result 
for the squared three-body amplitude in the soft limit as 
given in Eq. (A. 14) of Ref. [25] but with the Born contri- 
bution for 75 — > QQ being replaced by its spin-dependent, 
color-averaged counterpart, which reads in d = 4 — 2e di- 
mensions: 



1 



where 



2toqS 

t-\U-i 



- 1 



(7) 



(8) 



Here, t\ and u\ are the usual tree-level Mandelstam vari- 
ables. Contrary to the unpolarized case, AB jg receives no 
O(e) contributions [IB], cq is the electromagnetic charge 
of the heavy quark Q in units of the coupling e, i.e., 
e c = 2/3, and e 2 /47r = a. 

All 2—^3 processes exhibit singularities related to 
collincar splittings off the incoming photon (y — > +1) and 
parton b (y — > —1) which need to be factorized into the 
bare photon and proton PDFs at a scale fi / , respectively. 
As for the soft contributions discussed above, the kinemat- 
ics collapses to the much simpler case of 2 — > 2 scattering, 
such that compact analytical expressions, summarized by 
dAa^ in Eq. @, can be obtained: 



1 ~ \ ui 



(47r) £ - 2 r[l + e] 
1 



-l-e 



1-X 



2e 



4e 

log(l - x) 



-<ZPS 2 



1-x 



where 



Af\ c g -\x,O x ) = 32Tra s s(l - x) 

x A\M, 12 1 
Af£+\x,6 1 ) = 327rae 2 Q s(l-x) 



(9) 



(10) 



(11) 



X ^WqqV I AP qi {x), (12) 



Af^-\x,e 1 ) = m-Ka s s{l-x) 

AP gq {x). (13) 

Here, p\ and P2 denote the momenta of the photon and 
parton &, respectively. Both, the standard two-body phase- 
space dPS2 in © and the Born matrix elements squared 
in Eqs. (ITUl) - (fT3"]) arc to be evaluated with appropriate 
collinear kinematics as indicated by the shift of momenta 
P2 — > xp2, etc; further details can be found in Rcfs. 25, 

ESUTSJ. 

The resulting 1/e divergence in (JSJ assumes the form 
dictated by the factorization theorem, i.e., a convolution 
of (i-dimensional helicity-dependent LO splitting functions 
APij(x) and Born matrix elements Zi|M a b| 2 . By adding 



appropriate counter cross sections to ©, which to 0(aa 2 ) 
schematically read 



~2tt 



E 



dx 
x 



AV lJ (x 7 fif)dAa ( i l\xs) 



+ AV lb {x^ f )dAa { ^{xs) 



(14) 



where 



Vij{x,Hf) 



AP l3 {x){ 



In 47T - 



■\A] + A 9ij {x) 
A* 



(15) 

all collinear singularities can be consistently factorized 
into the scale evolution of the bare photon and proton 
PDFs, depending on whether they originate from collinear 
configurations involving the initial-state photon or parton 
b. The required Born cross sections and LO APij(x) are 
listed in the Appendix of Ref. [TH] , except for 



AP cn = C A [2x - 1 - 2e(l - x)] 



(16) 



Ca — 3 in (|T6| and, below, Cf = 4/3 are the usual QCD 
color factors 

The factorization scheme is fully specified by the choice 



of Agij in (fT5")) for which we take Ag q 



-4Cf(1 



to guarantee helicity conservation within the HVBM pre- 
scription for 75 in d dimensions [28], and Ag^ = other- 
wise. The Euler constant and ln47r, both, like the scale 
/i, artifacts of dimensional regularization, are subtracted 
along with the 1 je singularity. This defines the MS scheme 
in the polarized case, see, e.g., Ref. [25] . 

Finally, the last term in Eq. dAa^ , contains all 
the remaining, finite contributions, and the phase-space 
integration can be performed numerically in four dimen- 
sions. Monte Carlo integrations of ([3]) for different mea- 
surement functions S can be done in parallel by randomly 
generating a sufficiently large sample of final-state con- 
figurations characterized by xi, X2, x, y, 6i, and 02 to ac- 
count for the possibly large cancellations among the var- 
ious terms. Remnants of the regularization of soft and 
collinear regions of phase are the mathematical distribu- 
tions appearing in the different contributions to Eq. ([5]). 
Their proper definitions through a test function, proper- 
ties, and numerical treatment are discussed at length in 
Refs. [13123111] and need not be repeated here. 



3 Phenomenological Studies 

3.1 Preliminaries 

Based on our Monte Carlo code described above, we pre- 
sent a comprehensive study of charm quark photoproduc- 
tion in longitudinally polarized muon-deuterium collisions 
at a c.m.s. energy of V~S ~ 18 GeV relevant for the COM- 
PASS experiment at CERN. Unless specified otherwise, 
we will show results for single-inclusive cross sections, dif- 
ferential in the transverse momentum p® of the observed 
D meson. 



Johann Riedl et al.: Longitudinally Polarized Photoproduction of Heavy Flavors at Next-to- Leading Order of QCD 5 



The hadronization of the heavy quarks produced in the 
hard scattering is modeled using Eq. ((4]) with [30] 



D Q ^ H ®{z) =N Q z a Q(l-z) 



(17) 



The fragmentation function D®^ Hq is normalized such 
that its first moment is unity, i.e., Nq = (uq + l)(ag + 2). 
For the remaining free parameter in dTT)) we take a c = 5 
as the default value in all our calculations; see Table 4 in 
Ref . [31] . To estimate the uncertainties associated with the 
choice of a c , we will vary it in the range 3 < a c < 7 [31 . 
Like in the experimental analyses [15,22,23723]. we adopt 
a cut z > 0.2 throughout our phenomenological studies, 
where z is the fraction of the energy of the photon taken 
by the D meson in the laboratory frame. 

In the computation of the LO and NLO unpolarized 
cross sections we use the LO and NLO CTEQ6 parton 
densities [32] and values for the strong coupling a s , re- 
spectively. In the polarized case, we adopt the best fit 
from the comprehensive global analysis performed by the 
DSSV group [2j, which is the only set of helicity PDFs 
including constraints on Ag(x,/i 2 ) from pp data. The re- 
sulting PDFs are characterized by a small Ag(x,n 2 ) in 
the IE-range predominantly probed by RHIC experiments, 
J ° 2 r Ag{x, 10GeV 2 )dx = 0.005±g;i§|, with a node at x ~ 
O.i [2J. To study the sensitivity of charm photoproduction 
data to different Ag(x,/j,f), we also use the GRSV "stan- 
dard" set [33], which has a positive Ag that is larger in 
size than the one of DSSV. We note that the most recent 
RHIC data [3[9] tend to prefer a polarized gluon density 
somewhere in between the ones obtained in the DSSV and 
GRSV fits [HIE]- 

Equation ([3J consistently includes both direct and re- 
solved photon contributions to the spin-dependent photo- 
production cross section by defining Af l a as the convolu- 
tion of the polarized lepton-to-photon splitting function 
APji and the structure functions Af2 of a circularly po- 
larized photon, i.e., 



1 dy 

„ V 



AP 7l (y)Af2 



I) 



M/ . (18) 



Nothing is known experimentally about the Af2, which 
account for the hadronic structure of a photon in the re- 
solved case where a ^ 7 in (J3J. However, as will be demon- 
strated below, uncertainties associated with AJ2 turn out 
to be negligible for charmed meson photoproduction at 
COMPASS kinematics. Contrary to inclusive hadron pro- 
duction [18] , the resolved contribution turns out to be nu- 
merically small even when estimated with the "maximum" 
model for Af2 introduced in Ref. [33] and based on sat- 
urating the positivity limit \Af2(x, fi/)\ < f2{x,/j,f) at 
some low scale \x / using the unpolarized set of GRV pho- 
ton distributions [35] as reference. The direct part of the 
cross section J3J where a = 7 is obtained by setting 



Af2(x,n,) = 6(l-x) 



(19) 



in Eq. (fl8|). 

The collinear emission of a quasi-real photon with low 
virtuality Q and momentum fraction y off a muon with 



mass m M is given by the Weizsacker- Williams equivalent 
photon spectrum which reads in the polarized case [36] 



APM = Tn 



1 - (1 - y? ln Q 



2 

max 



(i-y) 



m 2 y 2 



2<y 2 



) 2 

- max 



i-y 

m 2 y 2 



(20) 



The upper limit Q ma x is determined by experimental con- 
ditions. For COMPASS we take Q^ ax = 0.5 GeV 2 and, in 
addition, restrict the fraction y of the muon's momentum 
taken by the quasi-real photon to the range 0.1 < y < 0.9. 

We use m c = 1.35 GeV as the value of the charm 
quark mass for all our results. For the factorization and 
renormalization scales in Eq. (|3]) we take /i/ = /i r = 
£,(Pt + rn 2 ) 1 ^ 2 with £ = 1 as the central value and where 
Pt denotes the transverse momentum of the charm quark. 
As is commonly done, we vary them simultaneously in the 
range l/2<£<2to estimate the residual scale depen- 
dence at NLO which represents the dominant source of 
theoretical uncertainty and can be taken as a rough mea- 
sure of yet unknown higher order corrections. 



3.2 Numerical Results 

We begin our detailed numerical studies with a discus- 
sion of the scale dependence of the unpolarized and po- 
larized photoproduction cross sections for D mesons in 
muon-deuterium collisions at COMPASS. Figure [2] shows 
the scale ambiguity from varying fj,f^ r = £,{p^ + m 2 ) 1 / 2 
simultaneously in the range 1/2 < £ < 2 for the dominant 
direct photon contribution (a = 7) only. The solid lines 
denote our default choice of scales, £ = 1. Since the pub- 
lished COMPASS data [TS] are divided into three bins in 
the energy Ed of the produced D meson, Ed < 30 GeV, 
30 < E D < 50 GeV, and E D > 50 GeV, we impose the cut 
Ed < 30 GeV in Fig. [5] Results for the other two bins in 
Ed are very similar and hence not shown. 

As can be seen, the theoretical uncertainties due to 
the choice of scale are quite sizable, both in the unpolar- 
ized and in the polarized case and can be even further 
inflated by varying [if and ii r independently as was done, 
for instance, in our study of heavy quark hadroproduc- 
tion 19J. Also, possible variations of m c , which we do not 
pursue here, would add to the theoretical error. By com- 
paring the solid and dashed lines in Fig. [2] one can infer 
the relevance of the resolved photon cross section (a ^ 7). 
As was argued in Sec. 13.11 its contribution is rather small 
compared to the direct photon cross section and potential 
uncertainties due to the unknown Af2 do not matter. 

Next, we study the role of the individual partonic sub- 
processes and their contribution to the photoproduction 
cross section in Eq. ([3J. Figure [3J displays the decompo- 
sition of the unpolarized direct photon cross section for 
£ = 1, as shown in Fig. [2J in terms of the 73 and jq pro- 
cesses. Again, we choose Ed < 30 GeV, and results for the 
other two bins in Ed are very similar. As expected, the 
photon-gluon fusion mechanism, which is already present 
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Fig. 2. Scale dependence of the single-inclusive transverse momentum spectrum of D mesons at NLO with Ed < 30 GeV for 
the dominant direct photon contribution in unpolarized (left-hand side) and polarized (right-hand side) muon-deuterium 
collisions at COMPASS kinematics; see text. Factorization and renormalization scales are varied simultaneously in the range 
(j, = £(py + m c) 1//2 > 1/2 < £ < 2; the solid line refers to £ = 1. In both panels, the dashed line shows the result for the sum of 
direct and resolved photon cross sections for £ = 1. 











do D / dp° [jib/GeV] : 










CTEQ 6 : 
m =1.35 GeV, E n < 30 GeV 

c ' D 






1 1 ■ 




'"■ NLO direct \ 

NLOyg J 

NLOyq ; 

L , LO direct [ 




J 











0.5 1 1.5 2 2.5 

PxfGeV] 



Fig. 3. Contributions of the individual 7<? and 75 initiated 
subprocesses to the direct photoproduction cross section at 
NLO for £ = 1 and E D < 30 GeV (solid line). Also shown 
is the LO result (dashed line). 



at LO, gives the main contribution to the cross section. 
The genuine NLO photon-quark channel yields a negative 
but small correction. The NLO result for the D meson 
production cross section is roughly a factor of two larger 
than the corresponding estimate at LO as can be gathered 
from comparing the dashed and solid lines in Fig. [3] 



A similar exercise in the polarized case is shown in 
Fig. 2] for two different choices of helicity PDFs. First, 
one notices that the results obtained with the DSSV and 
GRSV sets differ in sign and magnitude of the cross sec- 
tion, which is readily explained by the very different gluon 
densities in both sets. The positive definite Ag(x,/if) of 
GRSV leads to a similar decomposition into 73 and 79 
subprocesses as was observed in the unpolarized case in 
Fig. |31 Again, the cross section is strongly dominated by 
photon-gluon fusion, and the channel only yields a 
small correction. On the contrary, the oscillating Ag(x, ///) 
of the DSSV set of helicity PDFs leads to a negative 
dAo D . Since the DSSV gluon is much smaller in size than 
the one of GRSV, the genuine NLO photon-quark con- 
tribution, which is numerically very similar in both PDF 
sets, is more important and yields more than a quarter of 
the cross section at small p®- Another important obser- 
vation concerns the relevance of NLO corrections which 
appears to be very different for the DSSV and GRSV he- 
licity PDFs. This implies that higher order QCD effects 
do not cancel in the experimentally relevant double-spin 
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dAo D / dp° [LLb/GeV] 
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m c = 1.35 GeV,E D <30GeV 
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Fig. 4. As in Fig. [3] but now for the spin-dependent photoproduction cross section for two different choices of helicity PDFs: 
DSSV 2 (left-hand side) and GRSV "standard" [33J (right-hand side). 



asymmetry, 



which we shall discuss in more detail below. Clearly, for 
a reliable quantitative analysis of charm photoproduction 
in terms of polarized PDFs, preferably as part of a global 
QCD fit, it is indispensable to properly include NLO cor- 
rections. 

Figure [5] illustrates the uncertainties due to our choice 
of a c in the fragmentation of the charm quark into the 
observed D meson. As in Fig. [2] we show the direct pho- 
ton contribution for £ = 1, but here we do not impose 
any cuts on the D meson to allow for a better compa- 
rability with the cross sections on the charm quark level 
which are presented as well. Compared to the factoriza- 
tion scale ambiguities, the dependence of our results on 
a c is fairly moderate for both the unpolarized and the 
polarized cross sections. A similar observation was made 
in the case of hadroproduction [19]. Since charm quarks 
lose only very little of their momentum during hadroniza- 
tion, i.e., D®^ Hq (z) is peaked at fairly large values of z, 
the convolution ([4]) introduces only a rather small shift in 
the transverse momentum spectrum of the charm quarks. 
This can be inferred from the curves labeled d(A)<r / dp^ 
in Fig. [5] 

Another interesting question concerns the range of mo- 
mentum fractions x predominantly probed in the PDFs by 
the COMPASS data. Due to the dominance of the photon- 
gluon fusion process, charm photoproduction will mainly 
lead to a constraint on the gluon helicity distribution, 
which is the prime motivation for such measurements. The 
x distribution in LO and NLO, for both the unpolarized 



and polarized direct photon contribution to the cross sec- 
tion © is explored in Fig. [5] No cuts are imposed in this 
calculation and the D meson spectrum is integrated over 
the entire phase space. As has to be expected from the 
fact that the PDFs enter the cross section ^ through 
a convolution, a broad range of x values is sampled. It 
turns out, however, that the mean value of x, where the 
distribution d(A)a D jdx is peaked, is fairly independent 
not only of the order in perturbation theory, LO or NLO, 
but to a large extent also of the chosen set of polarized or 
unpolarized PDFs. We roughly estimate the average mo- 
mentum fraction to be (x) ~ 0.08 with an error of about 

+0.12 
-0.03- 

Our results differ from preliminary estimates of (x) by 
the COMPASS collaboration [22][2l[24], where NLO re- 
sults have been obtained based on some parton shower 
Monte Carlo to approximate the phase space for the NLO 
matrix elements of Ref. [16]. Significant differences be- 
tween (x) estimated in LO and NLO are found in this 
way. While their LO result for (x) agrees with our esti- 
mate of about 0.08, their preliminary NLO result is (x) = 
0.281" Since the details of the method are not yet pub- 
lished, it is not yet clear how these results can be com- 
pared to our full NLO calculation. We also note that once 
data on photoproduction processes are implemented in 
global QCD analyses of helicity PDFs, information on (x), 
though useful, is no longer required or relevant as the fits 
automatically impose the constraints from data for any 
given functional form assumed for the Af(x, /if). 

Next, we turn to the experimentally measured double- 
spin asymmetry, defined in Eq. (|2ip. which was analyzed 
at LO accuracy and under certain simplifying assump- 
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x 10 

0.01 I I 0.02 I 




Fig. 5. As in Fig. [2] but now for variations of the parameter a c in the fragmentation function (|17l) for the direct cross section 
computed with f = 1. No cuts on the D meson or charm quark are imposed; see text. Also shown are the results on the charm 
quark level, d(A)a /dp?. 



0.05 



0.03 



0.02 



0.01 



-0.01 



-0.02 




d(A)a u / dx [|j.b] 
NLO m c =1.35GeV 

da D /dx, CTEQ6 

dAc^/dx (x 50), DSSV 



dAa u /dx (x 50), GRSV 



0.05 



0.0.1 



0.01 



-0.01 



-0.02 




d(A)a u / dx [|xb] 

LO m c =1.35GeV 

da D /dx, CTEQ6 

dAo D /dx (x 50), DSSV 

dAa D /dx (x 50), GRSV 



Fig. 6. Distribution d(A)a D / dx in the momentum fraction x probed in the PDFs at NLO (left-hand side) and LO (right- 
hand side) accuracy for both the unpolarized and polarized direct photon contribution, integrated over phase-space and without 
imposing cuts. Note that the polarized results are scaled by a factor of 50 for better visibility. 

tions in terms of the mean gluon polarization Ag( (x), licity PDFs used throughout this paper, by computing 

g((x),fif) by the COMPASS collaboration [15]; see the the ratio of the cross sections shown in Figs. [2] and [3] for 
discussion below. Preliminary NLO estimates are also avail- Ed < 30GeV. To resolve the differences in All obtained 
able at NLO with the hybrid method outlined above [22J with the DSSV and GRSV PDFs, which mainly stem from 
[331121] . In Fig. [7] we show All for the two sets of he- 
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P? [GeV] 

Fig. 7. Upper panel: double spin asymmetry in LO and 
NLO for D meson photoproduction at COMPASS kinematics 
for the DSSV and GRSV sets of helicity PDFs; as in Figs. 
and [3] a cut Ed < 30 GeV is imposed; lower panel: the corre- 
sponding ratio of the NLO and LO results. The region where 
All changes sign and the ratio becomes meaningless is indi- 
cated by the shaded band. 
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Fig. 8. Predictions for the double-spin asymmetry for D me- 
son photoproduction at COMPASS at NLO accuracy for three 
different sets of helicity PDFs and in three bins of Ed com- 
pared to data. Note that in the bottom panel one of the data 
points is outside the range shown for All- 



Ag, an experimental precision of at least SAll ~ 0.02 
needs to be achieved. 

We also compare NLO and LO estimates of All for 
both sets of helicity PDFs in the lower panel of Fig. [7] and 
find rather different patterns depending on p®. At small 

, the NLO All is about a factor of two larger than 
the LO estimate for the DSSV set whereas a reduction 
by roughly the same amount is found with the GRSV 
PDFs. This illustrates that any approximations for the 
spin asymmetry, either to assume a cancellation of NLO 
corrections or a constant pattern independent of the choice 
of PDFs, are not justified and should not be used for an- 
alyzing data. Again, only a global analysis will lead to 
consistent results. Qualitatively very similar results have 
been obtained for other cuts on the energy Ed of the ob- 
served D meson, 30 < E D < 50 GeV and E D > 50 GeV. 

Finally, we compare our calculations at NLO accu- 
racy with the available results from the COMPASS ex- 
periment. Figure [5] shows the data in three bins of the 
energy Ed of the detected D meson as a function of its 
transverse momentum p^, imposing the cut z > 0.2. Note 
that instead of using the published data [15], we show 
new, preliminary results presented recently in Ref. 22 , 23 , 
[24] . A weighted average is performed to combine the re- 
sults for the three decay channels D — > Kn, D — > Ktttt, 
and D Rtttttt listed in [2"2H23U24] . We adopt the photon 



polarization dilution factors tabulated in Ref. 22,23,24 
to convert the data for All given in the photon-nucleon 
system to the double-spin asymmetries for muon-nucleon 
scattering computed with our Monte Carlo code. 

Clearly, present experimental uncertainties are too large 
to discriminate between different spin-dependent gluon 
densities. Apart from the DSSV and GRSV "standard" 
sets, we also use an extreme GRSV set, characterized by 
a very large and negative gluon density based on setting 
Ag(x,fio) = —g(x,fi ) at some low bound-state like in- 
put scale 33 for the evolution. Despite leading to a 
distinctively different All than our two default sets, all 
results are compatible with data within the experimental 
precision. 

To overcome the statistical limitations, fewer bins have 
to be used or all data need to be combined. COMPASS has 
performed such an analysis [T5J yielding the result for the 
gluon polarization quoted in Fig. [JJ which is compatible 
with similar extractions of Ag/ 1 g from hadron production 
data. These kind of analyses are LO estimates, assuming, 
in addition, that the convolutions of the PDFs with the 
partonic hard cross sections in Eq. ([3]) can be approxi- 
mated as 

All _ Ag g dA °ig M A 9 I dAa ig \ ^ 
g <g> dAa lg g \ da lg I 
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Here, the "analyzing power" ( d °~' 9 ) is evaluated at some 
average kinematics. As mentioned before, in a proper glo- 
bal analysis of helicity PDFs one does not need to resort 
to any of these approximations whose validity is difficult 
to estimate or justify. With our new Monte Carlo program 
it is now possible to perform such an analysis at NLO ac- 
curacy in the future. As we have demonstrated in some 
detail, NLO corrections are indispensable for a quantita- 
tive analysis. Most importantly, they do not cancel in the 
spin asymmetry as one might naively expect. 

Apart from the phenomenological applications presen- 
ted here, our code can be also used to estimate spin asym- 
metries for charm and bottom photoproduction and their 
impact on our knowledge of the spin structure of the nu- 
cleoli at higher c.m.s. energies intended for a first polar- 
ized lepton-ion collider. Details on the EIC project can 
be found in [10] . Such a machine will be indispensable to 
finally unravel the quark and gluon contributions to the 
nucleon spin [llj . 



4 Summary and Outlook 

We have presented a flexible parton-level Monte Carlo 
program to compute heavy flavor distributions at NLO 
accuracy in longitudinally polarized lepton-nucleon colli- 
sions in the photoproduction regime. For the first time, 
we consistently include both direct and resolved photon 
contributions and found the latter to be negligibly small 
for charm production at COMPASS kinematics. Experi- 
mental acceptance cuts, the hadronization of the produced 
heavy quark pair, and, if needed, their subsequent semi- 
leptonic decays can be included in phenomenological ap- 
plications. 

Heavy flavor photoproduction receives its importance 
for the field of spin physics from its expected strong sen- 
sitivity to the polarized gluon density which we confirm. 
In general, higher order corrections are found to be siz- 
able and strongly dependent on the chosen set of helicity 
PDFs. This is also true for the experimentally relevant 
spin asymmetry despite naive expectations that QCD cor- 
rections cancel in the ratio. Theoretical uncertainties due 
to the choice of the factorization scale are sizable even at 
next-to-leading order accuracy while ambiguities from the 
exact form of the charm quark fragmentation function are 
less important. 

The results obtained in this paper allow one to include 
data on charm photoproduction consistently into future 
global analyses of helicity PDFs. We have shown that data 
from COMPASS can lead in principle to a constraint on 
the polarized gluon density at a momentum fraction of 
about 0.1. Currently available data will have, however, 
very little impact on existing fits due to the size of the 
experimental uncertainties which are too large to discrim- 
inate between different gluon densities. 

In addition to the phenomenological studies performed 
in this paper, our code will be useful in assessing the 
physics impact of heavy quark photoproduction at a possi- 
ble future polarized lepton-ion collider like the EIC project. 



Note Added 

While completing our analysis, the preliminary COMPASS 
results [2~2l231l24] discussed in this paper have been pub- 
lished in [37| ■ O ur conclusions remain unchanged. 
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